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I investigate the vortex lattice structure of the Ginzburg Landau free energy for a two component 
order parameter in the weak-coupling clean-limit when the field is along the high symmetry axis in 
a tetragonal crystal. It is shown that the vortex lattice phase diagram as a function of the Ginzburg 
Landau free energy parameters includes phases with a hexagonal, centered rectangular, rectangular, 
and square unit cells. It is also shown that the square vortex lattice has the largest region of stability. 
The field distribution of the square vortex lattice near H C 2 is determined and the application of this 
model to Sr2Ru04 is discussed. 
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The oxide Sr2Ru04 has a structure similar to high T c materials and was discovered to be superconducting with a 
T c = 1.35 K by Maeno et, al in 1994 Q. It has been established that this superconductor is not a conventional s-wave 
superconductor: NQR measurements show no indication of a Hebel-Slichter peak in 1/T\T Q , T c is strongly suppressed 
by non-magnetic impurities ||, and tunneling experiments are inconsistent with s-wave pairing J9]. More recently 
there have been two experimental results that shed more light on the nature of the superconducting state. The /iSR 
experiments of Luke et. al. indicate that the superconducting state breaks time reversal symmetry which implies that 
the superconducting order parameter must have more than one component Q. Of the possible representations (REPS) 
of the point group, the two dimensional (2D) T 5u representation (REP) is the most likely state that exhibits 
this property. The order parameter in this case has two components (^7i , 772 ) that share the same rotation-inversion 
symmetry properties as (k x ,k y ) Q. The broken T state would then correspond to Theoretical 
arguments supporting a triplet pairing state have been given in Ref. Q. Given that this material may well be 
described by such an order parameter, it is of interest to explore further consequences of a two-component order 
parameter. In an earlier work it was shown that a consequence of the low temperature broken time reversal symmetry 
state is that the mean field vortex lattice phase diagram will exhibit two vortex lattice phases when the field is along a 
high symmetry direction in the basal plane [0 . Another important experimental development is the observation of a 
square vortex lattice in Si~2Ru04 by Forgan et. al ||. Within the context of the orbital dependent superconductivity 
model for Sr2Ru04 the orientation of the square vortex lattice relative to the underlying ionic lattice dictates which of 
the Ru orbitals exhibit superconductivity [^) . This work focuses on the magnetic field distribution and the structure 
of the vortex lattice for the field along the c-axis for this two-component model. 

The free energy for the representation of the tetragonal point group is given by Q 

/ = - a\ij\ 2 + (3x\7^/2 + faimrii - mV*i) 2 /2 + [3 3 \ m \ 2 \ m \ 2 + Kl (\D xVl \ 2 + \D yV2 \ 2 ) (1) 
+ K 2 {\D ym \ 2 + \D xm \ 2 ) + k 5 {\D zVi \ 2 + \D zV2 \ 2 ) 

+ K3[(D x rn)0vmT + h.c] + Ki[0 v Tn)0 w Tl2)* + h.C.] + h 2 /(8n). 

where a = ao(T — T c ), Dj — Vj — jr^Aj, h = V x A, and A is the vector potential. To simplify the analysis 
the Ginzburg Landau coefficients are determined within a weak-coupling approximation in the clean limit. The 
measurements of Mackenzie et. al. of T c as a function of impurity concentration show that the ratio of the mean 
free path to the zero-temperature coherence length is > 8 for T c > 1.3 K || indicating that the clean limit should 
be a reasonable approximation for Sr2Ru04. Without an experimental knowledge of the characteristic frequency of 
the boson responsible for the pairing (presumably ferromagnetic spin fluctuations) it cannot be determined that the 
weak-coupling limit is appropriate for Sr2Ru04. Note that the spin fluctuation theory of Mazin and Singh indicate 
that T c /Tp sa 10~ 2 where Tp is the characteristic paramagnon frequency [jll]. This estimate in conjunction with 
T c /Tf ~ 10 indicates that the weak-coupling approximation is reasonable for Sr2Ru04, but further experiments 
are required to ensure this. Taking for the T$ u REP the gap function described by the pseudo-spin-pairing gap matrix 
(note that this choice is not unique): A = iiriiVx/^v 2 )) 1 / 2 + rj2V y /{{v' 2 .)) 1 '' 2 ]a z ay where the brackets () denote an 
average over the Fermi surface and oi are the Pauli matrices, writing rj + = (rji + ir/2)/V2, = (^1 ~ ^Vv^, and 
rotating (D x ,D y ) by an angle 8 about the z-axis to obtain (D x ,Dy), the following dimcnsionless free energy is found 
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,/k — iA v , / is in units -B 2 /(47r), lengths are in units A = [H c 2 0\/(32e 2 Kian)] 1 S 2 , 
= 4^/(3 + i/),i/= ««*) - 3<t£t#)/«t£) + (« 2 z; 2 )), & is in units y/2B c = $ /(27rAC) 
(here B c has been chosen to represent the thermodynamic critical field), a — ao(T—T c ), £ = (Kw/Soe) 1 ' 2 , and k = A/£. 
The parameter j/ (note \v\ < 1) gives a measure of the square anisotropy of the Fermi surface. For a cylindrical Fermi 
surface v = and for a square Fermi surface = 1. It is easy to verify that in zero-field (rji, 772) oc (1, i) is the stable 
ground state for \v\ < 1. 

For the magnetic field along the c-axis the ground state is found by setting D z — 0. Writing II + = ^fn(iDx + 
Dy)/\/2H and II = ^fn{iDx — Dy)/y / 2H, minimizing the quadratic portion of Eq. || with respect to ?y+ and 
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where N = II + II_. The maximum value of the externally applied field H that allows a non-zero solution for (77+, 7?_) 
yields the upper critical field H C2 . At H = H C 2 the vector potential is that for a spatially uniform field A = (0, Hx, 0). 



Expanding , ?y_ ) in terms of the eigenstates of N (Landau levels) up to N 
matrix yields the result for en — n/ H C2 shown as a function of v in Fig. |l|. 
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FIG. 1. en = k/H C2 as a function of v. 

The form of the eigenstate of the H C2 solution is found to be r] + (r) — J2 n >o CL4n + 2( f >4n + 2 ( r ) anc ^ 7 ?-( r ) = 
E n >o a4,i ^ 4n ( r ) where </>n(r) are the harmonic oscillator wave functions (Landau levels). As is well known these 
wave functions have a large degeneracy, and the form of the vortex lattice is found by including the nonlinear terms 
of the Ginzburg Landau equations perturbatively to break this degeneracy. Following the procedure of Abrikosov, 
the average Gibbs free energy density is found to be (a derivation of this result for unconventional superconductors 
can be found in Rcf. E2| and for a mixed d and s-wave order parameter in Ref. fll3||) 
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fi is the fourth order homogeneous free energy, h s is the field (along the c-axis) induced by the supercurrent, and the 
over-bar denotes a spatial average. The form of the vortex lattice is found by minimizing (2k 2 — 1)(3a- To do this h s 
must be found. By minimizing the Ginzburg Landau free energy with respect to the vector potential the following 
relation is found for h s 
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Near H C2 the field h s is found by substituting the vector potential for the homogeneous field and the order parameter 



solution near H C2 into the right hand side of Eq. 
writing h s = J2 n m hn, m K<Kn yields (see Ref |l4 



Writing the left hand side of Eq. § as % — J2 n m Ui) 



h,l= J2n=lU+)n+l,n/Vn + T 

h,p = [Vt(j+)l-l,p - y/PU+)p-l,l]/(p- 1 ) 



and 
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The form of (3 a and k will therefore be determined by terms of the type n (r)0* i (r)0 p (r)0^(r). To evaluate such 
terms I make the assumption that the vortex lattice unit cell contains one flux quantum. The shape of the unit cell is 



kept arbitrary and the convention introduced by Saint- James et. al |15 to describe the unit cell is used. The lattice 
geometry is depicted in Fig. 0. 
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FIG. 2. The vortex lattice unit cell 



The lattice vectors are ai = a(l,0) and a 2 = fo(cosa,sina) with the single flux quantum constraint a6sina = 27r 
where a and b are in units Ih = 2tt-\/\£/H. This assumption allows the functions <p n to be written as fll2f 

Mr) = 2-/ 2 ^- 1 / Vr V2 E ^^(-^W'e-^'J 2 /^^ - y m ) (9) 

m 

where c n = e * 7r «(p+ 1 -"p) j y n — ( n _ \/2)^/2~tut, a = (b/a) sina, p = (b/a) cosa, and H n are the Hermite polynomials. 
To evaluate (3a and k it is useful to express the spatial averages in terms of a sum over the reciprocal lattice of the 
vortex lattice. The reciprocal lattice is given by G = iiki + where ki = - 2 ^ - (sin a, — cos a) k 2 = b 2 ^ a (1, 0). 
The general form of (3a becomes 
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where the coefficients a n ,-m,p,i and a n are determined by / 4 , h s , and the form of the eigenfunction near H£ , (f)h,i 2 — 
(a&sina) -1 J e~ lG r /(r), and mc denotes the unit cell. The following relation makes this formulation for determining 
(3a and k useful (found using the addition theorem for Hermite polynomials) 

2 -(m+p)/2 P ™ C rs H p J-z h l2 ]H m - s [z? , ] 
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The following relation is also useful 

{4> P 4>l)hM = [V2z W2 n|0 o | 2 )w 2 . (14) 

The relations 11-14 arc straightforward to implement numerically. 

In the analysis of the form of the vortex lattice the parameter v was considered to lowest order in perturbation 
theory (recall that \v\ < 1 so that v provides a natural expansion parameter) . The limit v = was considered by 
Zhitomirsky who analytically found the ground state eigenvector near H C2 p6| . The solution of the order parameter 
to first order in v is 

(V+,V-) = [0o + b^e-^fa, ~e m {efo + b 6 ve- m cj> e )] (15) 

where e = V3 - V2 « 0.31784, 6 4 = 2%/3e(10 + V6)/(36 + 16^6) w 0.18230, and 6 6 = V306 4 /(10 + V6) « 0.080203. 
Substituting this solution for the eigenstate near H C2 into the Eq. 6 yields the coefficients 
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Application of Eq. || yields 

h s = ^[(1 - 3/V2e + 2e 2 )|0 o | 2 + (2e 2 - e/^)|0i| 2 + e 2 |0 2 | 2 (16) 

+ ^{{VWebi/4 - V66 6 /4)(e l4 Vi05 + e" l4e ^0i) 
+ (x/30e6 4 /4 - e& 6 - V26 6 /4)(e j4 V^ 6 + e~ i4e <%<h) 

+ (x/3e/2 - M(e l4 Vl0o + e" i4 Vo04)]. (17) 

Using this expression for h Sl (2k 2 — l)/3 'a should be minimized with respect to 9, a, and p to find the form of the vortex 
lattice. It can be proven when v > (y < 0) (2k 2 — 1)/3a can be minimized for 9 = tt/4 (9 = 0). For v = 0, (2re 2 — 1)/?_4 
is independent of 0. This is to be expected since v = corresponds to a cylindrically symmetric Fermi surface. It is 
also found that k varies weakly (« 0.01) for the different vortex lattice structures studied in this article (such behavior 
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is also present for mixed d and s-wave order parameters 13 1). While this variation is small it determines the form of 
the vortex lattice in the region of k ~ 1/V2 and at small k the vortex lattice phase diagram becomes quite rich. 

The form of the vortex lattice found in the large k limit agrees with that found under more restrictive assumptions 
in Ref. S. In this limit the lattice structure depends upon v. The behavior of the vortex lattice as a function of v is 
similar to the behavior as a function of temperature found for borocarbide Jl7| and d-wave [ fL8|]T3"| superconductors. 
For v = a hexagonal lattice is found. As \v\ increases the lattice deforms continuously until \v\ = 0.0114. For 
< \v\ < 0.0114 the vortex lattice is a centered rectangular lattice as shown in Fig. 3 (which can be described by 
9 = 7r /4 — a and p — cos a and a — sin a where 7r/3 < a < 7r/2, a ~ 7r/3 corresponds to the hexagonal lattice 
and a = tt/2 to the square lattice). For \v\ > 0.0114 the vortex lattice is square. If v > 0.0114 the vortex lattice is 
rotated 7r/4 with respect to the underlying crystal lattice while for v < —0.0114 the vortex lattice is aligned with the 
underlying crystal lattice. 

As mentioned above when k becomes sufficiently small the vortex lattice phase diagram becomes richer. Fig. |3| 
shows the region of stability for the three vortex lattice states that were found to be stable. In addition to the two 
phases described above a third phase appears for small k. The vortex lattice for this phase has a rectangular unit 
cell and is described by p — and a — b/a. This phase is stable because k is smaller for this phase than for both 
the square and hexagonal lattices. For a cylindrically symmetric Fermi surface (v — 0) and k < 0.75 the hexagonal 
vortex lattice is no longer the stable structure. This arises because k is in a local maximum for the hexagonal lattice. 




FIG. 3. The vortex lattice phase diagram as a function of the Ginzburg-Landau ratio « and the square anisotropy parameter 
v. The phase diagram is the same for v < 0. For v > (y < 0) the square vortex lattice is rotated 7r/4 (0) with respect to the 
underlying crystal lattice. 

The Type I to Type II transition can also be determined and it is not given by k = 1/V2 but by H C2 = H c (which 
corresponds to n — e-^/3/2 up to corrections that are second order in v). For a conventional superconductor k = l/\/2 
and H C2 = H c are equivalent. Here it is found that the k for which k = 1/V2 is less than k = for all lattice 

structures studied. An analysis of the Gibbs energies indicates that the Meissner state is the stable phase for H near 
H C2 when n < e-y/3/2. 
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FIG. 4. Contour plots of the induced magnetic field h s for a square vortex lattice with v 
The contours from darkest to lightest correspond to nh s — —0.42 — > —0.042 in units 0.042. 



(top) and v — 0.2 (bottom). 



Clearly the square vortex lattice has the largest region of stability in Fig |[ To further investigate the square vortex 
lattice the spatial variation of the magnetic field as given by Eq. |l^ is determined. This is shown in Fig. [I| for v = 
and for v = 0.2. The induced field h s has (in addition to a global maximum and a global minimum) a local minimum 
and a saddle point. Fig. shows the field distribution for these two values of v as determined from 



P{h) 



Jd 2 rS[h~h(r)} 
Jd?~r ' 



(18) 



The peak in P(h) is due to the saddle point in the spatial dependence of h s . As v increases the saddle point value of 
h s moves away from the minimum value of h s resulting in a larger 'shoulder' in P[h) as v increases. 

Now I turn to an application of these results to Sr2Ru04. This requires a determination of the parameters v and 
k. The value of k as defined above is given by 



enH C2 
V2H C ' 



(19) 



where en is given in Fig. [j] and H c (H C2 ) correspond to the measured critical (upper critical) field [note that the above 
choice of k and £ also implies H C2 — $ /(2e,f/7r£ 2 )]. In principle the values for H C2 and H c given in Ref. ]2(| can be 
used to estimate k, however the sample used had a T c = 0.9 K which indicates that impurities cannot be neglected 
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(since T" lax w 1.5 K) so that the clean-limit approximation used here is not valid. Until measurements on cleaner 
samples become available these measurements will be used to estimate k. Using the values of H C2 and H c given in 
Ref. |^(| yields n = 1.2e# < 0.7. This implies that either the square or the orthorhombic vortex lattices will occur 
depending on the value of v. 
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FIG. 5. Field distribution P(h) near H C2 for a square vortex lattice with v = and v = 0.2 



To determine the value of v experimentally the anisotropy of the upper critical field in the basal plane can be used 



g C2 (a) = 1_ 
H C2 (a + b) 1 



(20) 



where H C2 (a) [H C2 (a+b)] is the upper critical field for the field along a [a+b] measured at T c . This anisotropy has not 
yet been determined so that a microscopic model must be used to estimate v. LDA band structure calculations p2 23 



reveal that the density of states near the Fermi surface is due mainly to the four Ru id electrons in the ti g orbitals. 
There is a strong hybridization of these orbitals with the O 2p orbitals giving rise to antibonding n* bands. The 
resulting bands have three quasi-2D Fermi surface sheets labeled a, f), and 7 (see Ref. p3[|). The a and (3 sheets consist 
of {xz, yz} Wannier functions and the 7 sheet of xy Wannicr functions. In general v is not given by a simple average 
over all the sheets of the Fermi surface. A knowledge of the pair scattering amplitude on each sheet and between the 
sheets is required to determine v flfj|Jll|| . It has been shown that the pair scattering amplitude between the 7 and 
the {a, 0} sheets of the Fermi surface is expected to be small relative to the intra-sheet pair scattering amplitudes 
if Sr2RuC>4 is not an isotropic s-wave superconductor [frOf . This forms the basis for the model of orbital dependent 
superconductivity in which, to account for the large residual density of states observed in the superconducting state, 
it has been proposed that either the xy or the {xz, yz} Wannier functions exhibit superconducting order. This model 
implies that that there are two possible values of v\ one for the 7 sheet {v xy ) and one for an average over the {a, p} 
sheets {v X z,yz)- Using the following tight binding dispersions 



e 7 
e a,f3 



e° — 2t-y (cos k x + cos k y ) — 4i 7 cos k x cos k y 



(21) 



a,0 



2t a> f3(cosk x + cos ky) ± \ it 2 a p(cosk x — cos k y ) 2 + l&t 2 a ^ sin k x sin k y 



and using the tight binding values of Ref. |Tl| for the 7 sheet (e°,t 7 ,i T ) 



(-0.4,0.4,0.12) and the values 



( e a,/3i ta,p, ia.p) = (—0.3,0.25,0.075) for the {a,/3} sheets yields v xy — —0.6 and v xz ^ yz = 0.6. These values of v 
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seem too large since they imply an anisotropy of a factor of 4 in H C2 (a)/H c , 2 (a + b). However v depends strongly 
upon the tight binding parameters used, for example taking (e° t 7 ,£ 7 ) = (—0.52,0.4,0.16) yields v xy — —0.08. The 
qualitative result that v xy < and v xz . yz > is robust. Physically v xy < because of the proximity of the 7 Fermi 
surface sheet to a Van Hove singularity and v xz . yz > due to the quasi ID nature of the {a, (3} surfaces p2| , p3[ . 
Assuming \v\ < 0.2 implies k sa 0.7 in which case there will be a square vortex lattice that is rotated ir/4 (0) with 
respect to the underlying crystal lattice if the pairing occurs on the {a, /3} (7) Fermi sheets. It is encouraging that 
Forgan et. al have observed a square vortex lattice in Sr2Ru04 Further experimental studies of the vortex lattice 
should provide useful information as to the nature of the superconducting phase. 

In conclusion a Ginzburg Landau theory for a two component order parameter representation of the tetragonal 
point group has been examined with a magnetic held applied along the c-axis. The vortex lattice phase diagram near 
H C2 was found to be rich with a square vortex lattice occupying most of the parameter space. The field distribution 
of the square vortex lattice was determined yielding predictions for ^SR measurements. Finally, the application of 
this model to S^RuO^ indicates that a square vortex lattice is expected to appear. The orientation of the square 
vortex lattice with respect to the underlying crystal lattice yields information as to which of the Ru 4<i orbitals are 
relevant to the superconducting state. 
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